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Abstract 

We compute holographically the dispersion relation for a hydrody- 
namic mode of fluctuation (the phonino) of the density of supersymmetry 
current in A/" = 4 SYM at strong coupling. The mode appears as a pole at 
low frequency and momentum in the correlator of supercurrents. It has 
a wave-like propagation, and we find its speed and coefficient of attenua- 
tion. 
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1 Introduction 

In [1, 2, 3] the author of this note and his collaborators have proposed to 
look at the finite temperature phase of A/" = 4 SYM; general arguments 
lead to the prediction that the dynamics of low-lying excitations would 
be well described by the hydrodynamic approximation, where the only 
relevant degrees of freedom are the densities of conserved charges. The 
conjecture was borne out by the computation of the correlators of con- 
served currents, that exhibited the expected analytic structure, namely a 
pole for small frequency and momentum corresponding to excitations with 
a difi^usive or a wave-like behavior. The coefficients entering the disper- 
sion relations of these excitations, called transport coefficients [4] , capture 
the low-energy dynamics of the finite-temperature phase. We computed 
the transport properties of the stress-energy tensor and the R-currents. 
The calculation was performed holographically, using the known corre- 
pondence between SYM in 4 dimensions and type IIB supergravity in the 
AdSs X S'^ background [5]. This means that the results were valid in the 
limit of infinite number of colors and infinite 't Hooft coupling. 
Our calculation has subsequently been extended in very many directions, 
mostly trying to make predictions in more realistic theories by deforming 
the original setup, e.g. breaking conformal invariance, supersymmetry, 
adding fiavors etc. (see [6] for a partial list of relevant references). 

In this paper we fill a gap in the analysis of the linearized hydrodynam- 
ics of A/" = 4 SYM: we consider the transport properties associated to the 
fiuctuations of the supersymmetry current. It was first pointed out in [7] 
that a supersymmetric thermal medium must have a fermionic collective 
excitation, that they called phonino; in that paper they discussed the case 
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of the Wess-Zumino model, and determined the properties of the phonino 
for moderately small temperatures, but could not discuss the ultrarela- 
tivistic case, as it was too complicated. As far as we know, there have not 
been other attempts at computing the phonino spectrum perturbatively. 
The outline of this paper is the following: in section 2 we discuss the 
prediction for the correlators of supercurrents following from the hydro- 
dynamical Ansatz; in section 3 we detail the holographic computation of 
the correlators, which is essentially the solution of the classical equations 
of motion for the gravitino in the dual geometry. In the final section 4 we 
summarize and briefly discuss our results. 

2 Supersymmetric hydrodynamics 
2.1 The phonino 

We follow the treatment of [8], to which we refer for a more extensive 
discussion. We consider a 4-dimensional supersymmetric field theory. The 
density of the supersymmetry charge should be included in the list of 
efi^ective degrees of freedom in the long-wavelength regime, by the same 
arguments that hold for any conserved charge and imply that it will relax 
with arbitrarily long time scale. In order to turn on a vev for the density of 
the supercurrent we may introduce a fermionic chemical potential fiajfta- 
In the equilibrium state, the supersymmetry algebra implies ^ 

{S?) = ~(5{TM-^r , 

{St)^^l3m,){j'l.f- (1) 

The linearized hydrodynamics deals with small fiuctuations out of equi- 
librium. We assume that the relevant degrees of freedom are contained 
in the density of charge 5*" = p" , and a constitutive relation is needed in 
order to express the other components of the current in terms of . The 
most general expression, linear in fields and derivatives, and assuming 
that no other vevs and chemical potentials are turned on, is the following 

The third term in this expression is fixed by the susy algebra, and V, e 
are the pressure and energy density in the equilibrium state. If the theory 
is conformal V = e/3; moreover 7' Si is related by susy to and it 
has to vanish; this implies that Ds — Da- At this order, then, there is 
only one transport coefficient that captures the dynamics of the fermionic 
fluctuations. The conservation of the current, dtp -f djS-' = 0, closes 
the system of equations and allows one to find the existence of waves 
propagating in the medium; they are solutions of the form 

pa{t) = 6""°°'" * {5^cos{kvst) + {ki ■ ■y'"y^)sin{kvst)^ p(0) (3) 

■^Sec appendix A for our conventions on the gamma matrices. 
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with a velocity Vs = ^, to be compared with the usual speed of sound 
Cs = The transport coefficient Ds determines the damping of this 

excitation, that has been called "phonino" . 

In presence of a supercharge density, other conserved currents can 
also be affected via non-linear terms in their constitutive relations. The 
R-current receives a contribution 

V 

R = ^Pl P 

(where the dots are for terms not depending on p) and the coefficient 
Ds enters in the R-current correlators; in [8] it is shown that such non- 
linear and non-derivative terms in the constitutive relations make the 
correlator decay more slowly in time, with a power law instead of the 
exponential decay characteristic of diffusion processes. The leading large- 
time behavior of R-current correlators was found to be 

where w = e + 'P, 7,, = i]/w, xr and Dr are the charge susceptibility and 
diffusion constant. All these quantities have been computed for Af — 4 
SYM, with the exception (up to now) of Ds- 

2.2 Supercurrent correlators 

The presence of the fluctuation mode (3) has to be reflected in the struc- 
ture of the correlators of operators with the same quantum numbers (this 
is the statement of the fluctuation-dissipation theorem). In particular we 
are interested in the retarded correlator of the supercurrents Si, defined 
as 

Gtf{!^ = I d'xe-^'--e{x'){{St{x), 5f (0)}) . (5) 

The correlator (5) is constrained by the Ward identities: 

k'G^j = , -y'Gi, = G,j 7^ = . (6) 

The first equation is simply the conservation of the current, while the sec- 
ond follows from superconformal invariance as already mentioned. The 
correlator can be built out of projectors, that automatically implement the 
identities. First we analyse the case of zero temperature, where the corre- 
lator has to respect the 4-dimensional Lorentz invariance. The projector 
on the transverse gamma-traceless part of a spinor-vector is 

P.' (7- - ^) 7^ - ^ (4fc - 7.^ ■ (7) 

The correlator can then be written as Gij = Pi^ MkiP'' j- The matrix M 
can be expanded on a basis of gamma matrices, and the only form allowed 
by the symmetries is 

Mu = Aik") Vr^ki ■ 

The zero-temperature correlator is then completely fixed up to a scalar 
function of momenta, A(fc^). 
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At finite temperature, only spatial rotational invariance is preserved, 
and the correlator can depend, apart from k, also on another vector that 
is the velocity of the fluid; since it is at rest, this is it' — (1,0,0,0). Let 
us write fc* = (oj, q). 

We can form another projector with the properties of (7), that is also 
transverse to u^: 

Pii = Pu=0 

Pju = 5,u^\ (7. - ^)7fe - ^ (39. - lj<i)qk (8) 

and define = P - so that P^P^ = 0. The correlator can have 
several possible structures: , P"^ P'^ , and P-^P^ + P^P^. We are 

interested in the fluctuations of the density, which are contained in the 
longitudinal part, so we can write 

Gij = PikMkmPmj (9) 

and now the form of M allowed by the symmetries is 

Ckm ^ a{Vkm + bUkUm) , (10) 
C'km = a'iVkm + b'ukUm) ■ 

Not all the functions a,b,a',b' are independent, it turns out that the 
following combinations vanish inside the projectors: 

,,Ll) , , 

-¥^UkUm + y (Vkm + TT^UkUm ) ■ 

q-' 2q^ 

Thus for generic values of the momenta we can fix two of the parameters. 
We haven't been able to find a particularly nice canonical choice of struc- 
tures. The correlator is then determined, in this channel, by 2 functions of 
q, ui. The hydrodynamical prediction is that these functions should have 
a pole, in the complex lu plane, that goes to the origin as q ^ 0. 



3 Holographic derivation 
3.1 Zero temperature 

The two-point correlator of the supersymmetry current in A/" = 4 SYM at 
zero temperature has been computed in [9, 10]. We recall in this section 
the main points of the derivation, to be used for reference in the next 
sections. 

The operator dual to the susy current is the gravitino, a Rarita- 
Schwinger field propagating in the dual AdS geometry. We follow here 
the presentation of [11]. The bulk action is 

s= f d''+^xy^(*^r'*''''P),*p-m*^r'*''*.) , (11) 
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where r** = ea7"- For the time being we can leave d and m unspeci- 
fied. The action is evaluated on the (Euclidean) AdSd+i background with 
metric 

ds^ = \{dxl + S'^dxidxj) . (12) 

The study of the solutions of the equations of motion in the neighborhood 
of the boundary xq = shows that components of diflFerent chirality have 
different asymptotic behavior. In terms of the spinor xpa = £^''5'^ one has 

ipt{xo,Xi) = (a::o)''''^~'"(^i(x) + . . . , 
V-(xo,x) = (xo)'/'+'"x.(x) + ... (13) 

The solution is required to be regular in the bulk. It is possible to 
show that this imposes a relation between ipi and Xij so that only one of 
them can be assigned as a boundary condition. Without loss of generality 
one can take m > 0, then the appropriate boundary spinor can be shown 
to be ifi. Moreover, it has to satisfy a constrain 7' ifii — 0. The complete 
solution can then be expressed in terms of the boundary data : 



ipi{xo,k) = Xq^^ 



2k 1 



ixoh h {—;—)Km-l/2 



A:2 {2m + d~l)K^+i/2-2xokK^+s/2 
Here the Bessel functions have argument K,^{kxo). A similar expression 
holds for the conjugate field tp. 

We recall here one of the basic tenets of the AdS/CFT correspondence: 
the on-shell action, with given asymptotic conditions for the fields, is 
the generating functional for the correlators of CFT operators. The bulk 
action vanishes on-shell for fermionic fields, therefore the only contribution 
comes from the boundary action [12] 

Sidy = / d'':iiVhh"^,<lij . (15) 



The action has to be 
regularized by evaluating the integral on a surface xq = e, and taking 
the limit e — > 0. In the limit, the expression contains divergencies that 
in principle should be subtracted by introducing suitable counterterms, 
however one can see that the divergent pieces are analytic in k, therefore 
they only affect the correlators by contact terms. Finally, the regularized 
action can be written as 

Stdy = J dkS'' {^r{-^k)xAk)+Xr{-k)^Ak)) (16) 

the boundary value tp is the source for the field theory operator, and the 
appearance of a non-trivial correlator is due to the nonlocal relation that 
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expresses x in terms of ip, as implicitly given in (14). The result is, up to 
an overall factor, 



^ . 2(2m + 1) krks 

^ ~ (rf+2m + l)~fc2 



G., = ni (,?.. - , ^) n- , (17) 



where 11^ = (5[ — ^717*^ is the projector on the gamma-traceless part, since 
the source is gamma-traceless. The correlator is transverse, k^dj = 0, 
for m = that is when the gravitino is massless, as it should be. In 
the case of interest, d = 4 and m — 3/2, in units of the AdS radius. 

3.2 Finite temperature 

3.2.1 The gravitino in the near-extremal 3-brane back- 
ground 

In order to compute the correlator at finite-temperature, we need to study 
the Rarita-Schwinger equation in the background of a non-extremal black 
3-brane. In this case we are not able to solve the equations for generic 
mass and dimension, so we work directly in d = 4 and m = 3/2. The 
near-horizon geometry has the following metric : 



^ T R (jfy^^f _^ ^^2 ^ ^ 2 ^ ^^2N ^ R ,^2 ^^g^ 

u ^ ' 4/(?i)u2 

where /(«) — \ — 11} . 

The radial coordinate u £ [0, 1]; the boundary is at u = and there 
is a horizon at u = 1. Near the boundary, u is related to a;o of the 
previous section by u = x%. The covariant derivative on a spinor, in this 
background, is 

-Dp = + -C-jJlSab , (19) 



A = 9. + ic.f^5(0, (20) 

where i = t,x,y,z and (i) = 1,2,3,4. The non- vanishing components of 
the spin connection are 

51 7VT{1+U^) 

— 



52 53 54 m 
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The metric ( 18 ) satisfies i?^^ — Ag^,y with A — —A/R^ 
The Rarita-Schwinger equation is: 



(21) 
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The equation is invariant under a gauge transformation parametrised by 
an arbitrary spinor e: 

S-ii^ ^ Df,e+^mr^€. (22) 

We show in appendix B that the RS equations imply, in d space-time 
dimensions, 

For a generic value of the mass, (23) implies F ■ = 0, but for d = 5, 
A — and mR = 3/2, the coefficient in the brackets vanishes, so the 

constraint is not implied by the equations of motion. However, it can still 
be imposed as a gauge-fixing condition. We will work in this gauge. The 
relation (42) then implies that D • = 0. In holographic computations it 
would be perhaps more natural to impose an "axial" gauge = 0, but 
we were not able to solve the equations of motion in that gauge; at the 
level of two-point functions it does not make a difference, but it could be 
an issue in the computation of higher-order correlators. 
The equations of motion in this gauge can be written as 

F''(i:)^*„-£)„*^)+m*, =0 (24) 

Explicitly, and again redefining the field as = e^*!',!, eqs. (24) read 

(^' = a/9„v): 

21*2 — 35, It 1 , rriR , „ 
2uf f 2uJf 



— 2 5 , 1 +u'^ I , mR , „ 



+ ^7V.-^7.^. + ^^.^0 (27) 
j,k = x,y,z 

Even though we have fixed a gauge, there is still a residual gauge symmetry 
parametrised by e such that Pe H — = 0, or 

Notice that this equation determines only the radial profile of e, leaving 
arbitrary dependence on the spacetime coordinates. 
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3.2.2 Solution of the equations of motion 



The equations (25,26,27) are a complicated system that cannot be solved 
exactly. We want to solve them in perturbation theory in the momentum. 
We consider a fluctuation with a definite momentum fc'' — 2nT (ro, 0, 0, q). 

It is useful to decompose the fluctuations of the fields according to 
their spin under the unbroken transverse 0(2) (rotations in the xy plane). 
There is one component of spin 3/2, given by 77 = 7^^2 — 7'^ '03, and four 
components of spin 1/2, given by ip\, Tp4, Tpb and (f> = 7^^2 + ^^4>3. Com- 
ponents of different spin are decoupled, and we are interested in the sector 
of spin 1/2, since it contains the time component of the gravitino, which 
is dual to the supercharge density on the boundary, and we expect the 
correlators of the charge density to exhibit the hydrodynamic behavior. 
In this sector we can then set rj — 0. Moreover, we can use the matrix 
7^'^ as a chirality matrix; it commutes with the evolution, as one can see 
by inspection of the e.o.m., therefore we can choose to study the part of 
definite 7^'^ chirality; the number of components is then reduced by half. 

A further simplification comes from the fact that the gauge condi- 
tion 'j^jps + 'y'i'i = 0, together with the equations of motion, yields the 
following condition (compare with (2.17) of [9] ): 



3 \ i . /2u^/^ 1 Aim 



" " 4 ^^^^ 

where f = —ixx)l\/~f^^ -f 197''. 

We can use the relation (29) to solve for 7 ■ ^ = 7*1/)^ in terms of 
i\>\ and i\>4. Then we derive a system of equations for the two latter 
fields. Again we consider the 7^'^-chiral part. Finally, we have a system 
of 4 scalar equations for the remaining components that we write as a 
4- vector X — (tpf (here the superscript refers to the usual 

chirality in 4d). 

The resulting system has regular singular points at the boundary and 
at the horizon. The local solutions at the horizon are 



3 t» 

(1-^) 3 (i, 1,0,0) 

^i-uyi-'^ (0,0, i,i) 

(l-«)~'+ ~ (-i, 1,0,0) 
(l_u)-4+>- (0,0, -i,l) 

The first two eigenvectors corrispond to the incoming-wave boundary con- 
dition. As the previous finite-temperature computations have shown [2], 
one has to impose incoming-wave boundary condition at the horizon in 
order to recover the retarded correlator. The form of the eigenvalues 
shows that we can eliminate the outgoing wave solutions by imposing the 
boundary condition (1 -I- i75)'0i,4 — 0, or = iip~ . 

The solution at order zero with these boundary conditions is 
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3/4 



/ in«/*(2-y7)(l + \/7)-^/*«o \ 

V ^u^/^V7(l + ^/7)"'^'(^^/7"o + A.) / 
The solution depends on two parameters, ao and /3o- It can be checked 
that under the action of the residual gauge transformation (35), 

SctQ — ja+, Sf3o — — a+. 

The combination ao + i/3o is gauge invariant. We could use this freedom 
to fix one of the two parameters to zero, however we find it convenient 
not to do so and to keep the gauge invariance. At the next order, we 
have an explicit solution ^(i), whose form is very complicated, and is 
written in appendix C. The important point is that in solving the equa- 
tions in perturbation theory, at every stage there are 4 new integration 
constants to be determined. Again two of them can be fixed by imposing 
the incoming-wave condition at the horizon. A natural choice would be to 
require that each component has a prescribed behavior near the horizon, 
namely 

^^/-'>»/2(l + F(it)) 

with F a function of order 1 in momenta, and such that F{u = 1) = 0. It 
turns out that this requirement is incompatible with the gauge symmetry. 
However, the two remaining integration constants can be reabsorbed into 
a redefinition of qq and /?o to first order in momenta. It is not necessary 
then to keep track of them explicitly if we do not assume that a, [3 are ho- 
mogeneous functions of the momenta. The following non-obvious change 
of variables is convenient 

_ gp -I- iPo (q - 6ro)ao - iq/3o , , 

" ^ 2 ' 2 ■ ^ ' 

The boundary value of the spinors is given in terms of these parameters 
by 

% )~2(2-)^/^ ( (,V2+fq + 2a,(5-X+^))a+^/3 

(32) 

where L — log (l + \/2) • Solving for a, f3 in terms of (p\ , ip4 yields 



(3q-K>)yi-(q-3to)y4 , 
2~^/* I Pi^) I (33) 



where 



P(q, t») = -\/2j(q - 3t») - y q^ + 6to^ + 2\/2Lto(q - 3to) . 



This polynomial has a zero for small momenta at 

q 4\/2i 2 4n 
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This pole corresponds to the propagation of a mode with speed v — ^ 
and an attenuation given by the imaginary part of the dispersion relation. 
The value we find corresponds to a diffusion constant 

27VTD, = ^ (34) 

This is our main result. 

The reader may have noticed that while a was supposed to be a gauge- 
invariant quantity, its solution in terms of the boundary value is not easily 
recognized as such. However, just as in the zero-temperature case, the 
boundary value of the spinor is traceless, and the residual gauge trans- 
formation preserves this property. From (28) we find that the asymptotic 
boundary behavior of the gauge parameter is 

-1/4 1 •,/ 11/4 

e+ = a+u — —iijia-U , 

e- = a-u^/"*- iz^a+u'/*, ^ = -n>7^ + C)7"* (35) 
The transformation (35) acts on the boundary values as 

Sipl oc (q — 3t»)a+ (36) 
5(f4, oc (3q — vo)a+ 

and the combination appearing in a is invariant under this transformation. 

In order to compute the correlator we need to find the negative- 
chirality component of the fields at the boundary. The asymptotic be- 
havior of the negative chirality part is 

3/4- , 7/4 

and since in the boundary action the determinant of the induced metric 
gives a factor of and (p ~ m^^*, one can see that x contributes to the 

divergent part of the action, that has to be subtracted in the regulariza- 
tion. The subleading term x is the one that contributes to the finite part 
of the action. It is given by 

The other components X2, X3 s-i'g determined by the algebraic constraints 
= and 7^x2 - 7^X3 = 0. 
We can write the solution as a sum of two terms that depend on a 
and P respectively. Let us consider the a-dependent part: it depends 
on tpi only through the gauge-invariant quantity a, and it is transverse: 
fc'xi = 0; since the source is gamma-traceless, it is clear that inserting 
this term in the action (16) yields a correlator with the right properties, 
namely a conserved and gamma-traceless one. We can also put the result 
in a covariant form. We need to repeat the analysis for the components 
of the fields with the opposite 7^'' chirality. The equations are the same 
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after the exchange 1/14 — » —tp4,<\ — q. At the end we find an expression 
for the correlator that matches the general form (9,10), with 

a = 0, b' = I (38) 

81 i<f'-w'r 

2V2 (<l^ - 9tB2)2 • 

The second term in (37), depending on /3, is more troubling since 
it is gauge-dependent. However it does not contain the hydrodynamic 
pole, and under a gauge transformation the denominator is cancelled; this 
means that the contribution of X/3 generates a correlator that satisfies 
the Ward identities up to contact terms. In general we would expect 
such contact terms, on the basis of the following argument: the susy 
transformations in the bulk do not act only on the gravitino, but also 
on the metric (and on the gauge field, although this is not relevant for 
us). The susy variation of the vielbein is ^e^ — e'y""i! ^. When we try 
to prove the Ward identities, we should also include the boundary terms 
that depend on the metric; the boundary action is then : Sbdy = Sq + S2, 
where So contains the Gibbons-Hawking term and the boundary volume, 
and 5*2 is the term quadratic in the fermions that we have been using so 
far. Then under a susy transformation we have 

5S , _ . SS a T 

Taking the derivative with respect to 4', and then setting 4* to zero, we 
find the modified Ward identity 

{d,S^ S') + y ((7 • S)S') + (7O {r') = (39) 

However the precise form of the contact terms we find does not match the 
expectations, since /? only depends on ipi whereas the expected form, as 
appears in the formula above, looks like {T''-'}ynpj and so it depends also 
on (^4. It is possible that one should add some additional counterterms to 
the regularized action in order to restore the correct Ward identities [14], 
but we will not attempt to find them in the present paper. Nevertheless, 
we can still extract some information from the /3-terms. 

Given the structure of the boundary action, one sees that x as a func- 
tion of ip can be directly interpreted as the expectation value of the su- 
percurrent in the presence of sources. In the hydrodynamic regime, they 
should then satisfy the constitutive relations (2). It is possible to check 
this on our solution, in an expansion in momenta. Notice first that Xl3 is 
of order -1, whereas a and Xa start from order 0. Then the constituent 
relation at the lowest order, that is -1 in our case, reads 

and it is satisfied by the second term of (37). At the next order, from (2) 
we derive 

S{o) = g'S'Jo) +iDs (27rrq)Sf_i) 
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and one can check that this relation is satisfied, again only up to contact 
terms, with Ds given by the expression in (34). This provides an alterna- 
tive and more direct derivation of the value of the diffusion constant, more 
similar in spirit to the approach of [15], though of course at the linearized 
level the two methods are equivalent. 

4 Summary and discussion 

We have shown that the holographic description of a supersymmetric 
gauge theory plasma allows to exhibit the presence of a collective fermionic 
excitations, the phonino, whose existence can be predicted on general 
grounds; this excitation propagates like a sound wave, with a characteris- 
tic dispersion relation; the holographic method allows us to compute the 
lowest order terms in the dispersion relation, related to the speed and the 
attenuation. The results are valid in the limit of strong coupling. We find 
a finite value for the supercharge diffusion constant in this limit. The same 
was true for the other transport coefficients, like the shear viscosity and 
the R-charge diffusion constants, whose perturbative estimates indicated 
a vanishing value at infinite coupling. To our knowledge a perturbative 
calculation of Ds in a gauge theory has not been done. In this paper we 
have studied only the longitudinal part of the correlators, that involves the 
density fiuctuations. The transverse part could also be considered, and 
although it should not contain propagating modes, it should be possible 
to extract the value of the diffusion constant via a Kubo formula (see the 
footnote 32 in [8]). This would provide yet another confirmation of our 
result. 

It is tempting to speculate that Ds may have a universality property 
analogous to the shear viscosity, namely it has the same value in all su- 
persymmetric strongly coupled theories with a dual gravity description. 
It should not be too difficult to prove this using the methods of [15] or in 
the membrane paradigm (see [16]). 

Acknowledgments 

I owe a great debt to Andrei Starinets, who collaborated in the early stages 
of the project. I thank C. Bachas, J. Troost and K. Zarembo for useful 
discussions, L. Yaffe for pointing out an imprecision, and the Perimeter 
Institute where part of the work was done. The work was supported 
partly by the EU under the contract MRTN-CT-2004-005104 and PIEF- 
GA-2008-221026. 

A Gamma matrices 

We will use the following explicit representation of the five- dimensional 
flat gamma-matrices as 2x2-block matrices : 



13 



- ^ ( ; -I ) ■ («> 

where cr* are the usual PauU matrices. They satisfy {7^,7^} = 2riij with 
rjij — diag{—l, 1, 1, 1, 1). Notice that we call 1 the time component. 
In the text we often write a gamma matrix acting on a chiral spinor, it is 
then understood that the matrix should be replaced by the corresponding 
block with the right chirality; {■y^^)a = {o'^)ap'^^ ■ Chiral indices can be 
raised and lowered with e„^, e^^. 

B Gamma-tracelessness 

We follow the computation in [13]: Using the identity V^p" = rTT"^ - 
gMPp<^ - pP-^r^ + g^T", and that [D^J^] = 0, we find 

T^P{T ■ *) -_D''(r- *) -r"i3- * + ^*'' -mr"'"*^ =0. (41) 

Contracting the last equation with F^j wc find the relation 

PiV ^— i mr ■ * . (42) 

This can be plugged back in (41) to give, in d = 5, 

r (i3,*^ - D^*.) + ^ rT^*, + ^ r^r"*, = o . (43) 

Applying _Dp to (41) we have 

PP{T ■ *) - D^{r ■^)~pD-'f + D^P^^ - mr^"!)^*^ = . (44) 

Using PP — = ir^'^[_Dp, D^] one can rewrite this as 

ir'"'[D^, D,]V ■ * + r"[Li^, - m{pr • "I- - D • *) = , (45) 

Using once again (42) and that the curvature tensor acts on a spinor as 

l^ai,. 1 



we find 



\8 4 J d-2 

(46) 

After some gamma-matrix algebra, one can see that the expression in- 
volves in fact only the Ricci tensor 

R^y — Q R^pau — Af/^/y . 

Finally we obtain (23). 
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C Solution of the e.o.m. to first order 



Xi = 

{-/3/4(i + ^)l/4' 

(Ci- 



,9/4 



-2 + Vf) 



1 . /' 24yi+^(8q(-l+yj)+3(2+^)a>)ao ^ 



2 (^18t»Log[/] + 32x/3qLog [2 - 7/] - 72t»Log [l + V^ + Vf] - 
18^/2qLog 



27y2roLog 
64V3qLog 

64y3qLog , 

64qVl+v7(2ao-i)3o) 



V3 + vT+VJ\) 



ceo 



36qVl+yj(ao-i,3o) 



-1+V7 



+ 



9V2qLog 

^,,7/4 ^ 

^3/4' 



(C2 + 

1 /' 216Vl-V7t»Q 



1+^ 



(5ao - i/9o) ■ 
(2 + 77) 

96in(2q+3tX>)ao 



8g(ao+»/3o) 



+ 



Vi+v7(2+v7) 

72 v^( - i+ V^) lOLog [l - 0W7| «0 
2i+v^ 

18i%/2(2q + 3t») (Log [l + ^/7] - 2Log [^/2 - V^-Vf]) ao+ 
72it»Log [l + VT^TT] «o + M(£^o±^M_ 



36q^/l-v7(^Qo+;^o) 
i+v7 



\AW7(2+^) 

9y2q (Log [1 + ^ - 2Log [V2 - ^1^77] ) (5iao + A)) ) 

72/1/4 

(72Vl + V7(C3-iCiy7) + 
(9x/2\/l + VT{q - 6n))Log [x/2 - V^ + Vf] - 
3 (92q + 24t» + 3\/2\/l + ^//(g - 6t»)Log + VlTVf] ) 4 
^ (-180q + 72t» - 36\/rT77«»Log[/]- 
64x/3qVl + y/Log [2 - ^/J] 



9y2qVl + y/Log x/TT77 
54^/2/^T77tt>Log [^2 - x/TTTj 
64y3q/rT77Log [^3- x/TTTl] + 
144 V 1 + \/7t»Log [l + V^ + Vf 
9\/2qv/rT77Log [\/2 + ^1 + ^/7 
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54^/2y^T77t»Log [^/2 + \/l+ v^] + 

3i (^-20q + 12q\/7 - 48t» + 3^29 (-1 + ^/7) \/l + ^ 
Log [72 - Vl + yj] - 24^1 + V/wLog [-1 + \/l + v/j] + 
24^1 + \/7n'Log [l + v/TTTJ] + 
3V2qVl + %/7Log [^+ v/TT77 - 
3%/2qy7(lT77)"Log [v^ + ^/TTVT] ) A)) , 

1 -..7/4 



18/1/4 (l+v7)''' ™ ^ 
(I8C4 + 

i (18C2VJ - 18^^/7t»Log [1 - v/l-\/7] ao+ 
ISiVJroLog [1 + \/]~77] Qo- 

9(-i+V2)v7Log[l+v7]((q-6a))c,o-iq/3o) 

9 (-i + ^Af Log [^/2 - \/r^77] 

|o)_- - y7((31q - 126tt))Qo - 23iq/3o) + 



{{q - 6rD)ao 



4 (-i (-25q + V7(19<l - 36t») + 27t») ao+ 
(g (2-8/7) -9t»)/3oK 

18 ( -i^/2(q- 6t»)Log 



roLog 



i-Vi-v7 



- \/2qLog 



y2+\/i~yj 
\/i+v7 



A) 
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